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SECTION  1 


INTRODUCTION 


It  has  been  long  recognized  that  a  number 
of  ionospheric  processes  belong  to  the 
same  class  of  interchange  instabilities  as 
those  that  impact  nuclear  phenomenology. 
These  include  the  Rayleigh-Taylor  instabili¬ 
ty  that  operates  at  the  magnetic  equator,  the 

mid-latitude  fi  X  B  instability  that  is  driv¬ 
en  by  neutral  winds,  and  the  high-latitude 

E  X  B  instability  that  is  driven  by  magne- 
tospheric  electric  fields.  Large  barium 
clouds  released  into  the  ionosphere  also  are 
unstable  to  these  processes,  and  a  number 
of  experiments  have  been  performed  in 
space  to  study  them.  The  local  linear  theory 
for  these  instabilities  has  been  worked  out 
by  many  investigators  (see  review  by  Fejer 
and  Kelley.  1980.  and  references  therein). 

The  concept  of  “freezing"  associated  with 
plasma  irregularities  has  become  a  focus 
for  the  DNA  community.  Evidence  of 
freezing— that  is.  the  appearance  and 
persistence  of  F-region  plasma  irregulari¬ 
ties  at  a  preferred  scale  size— was  first 
noted  in  photographs  of  striations  produced 
by  the  CHECKMATE  nuclear  detonation 
|W.  G.  Chesmit .  personal  communication!. 
Since  then,  this  phenomenon  has  been 
observed  in  various  plasmas  at  F-region 
altitudes,  notably  barium  ion  clouds,  and 
the  nighttime  equatorial  and  polar  iono¬ 
spheres.  Freezing  characteristics,  and  the 
processes  that  produce  them,  are  of  interest 
because  they  appear  to  be  intimately 
associated  with  the  convective  interchange 
instabilities  and  are  expected  to  be  promi¬ 
nent  in  the  nuclear  environment.  Despite 
their  acknowledged  existence,  very  little  is 
known  experimentally  about  the  freezing 
characteristics.  The  freezing  scale  has  been 
observed  to  have  a  nominal  value  of  several 


hundred  meters  and  was,  at  first,  thought 
to  vary  little,  regardless  of  plasma  condi¬ 
tions.  This  nominal  freezing  scale  has  been 
determined  from  optical  measurements  of 
CHECKMATE,  of  barium  ion-cloud  stri¬ 
ations,  and,  as  will  be  presented  in  detail 
below,  from  spectral  analysis  of  in  situ 
plasma  density  fluctuations  and  phase 
scintillation  data. 

The  need  for  further  experimental  charac¬ 
terization  of  the  freezing  process  is  crucial 
because  relatively  untested  theories  for  its 
explanation  are  being  implemented  in 
system  codes  (e.g.,  SCENARIO,  MICE,  and 
MELT).  Moreover,  it  is  not  obvious  that  a 
single  mechanism  dominates  over  all 
others,  or  whether  two  or  more  processes 
contribute  to  the  formation  of  the  freezing 
scale.  It  is,  therefore,  particularly  important 
to  experimentally  test  candidate  theories. 
To  do  so,  we  have  examined  changes  in  the 
freezing  scale  under  varying  plasma  condi¬ 
tions;  we  examined  different  strengths  of 
instability  driver,  looked  at  contiguous,  but 
oppositely  directed  plasma  density  gradi¬ 
ents,  and  measured  changes  in  freezing 
scale  as  a  function  of  altitude  (collisional- 

ity) 

Implicit  in  the  following  analysis  is  the 
assumption  that  the  freezing  scale  (origi¬ 
nally  coined  to  describe  photographic  data) 
is  uniquely  associated  with  the  break  point 
between  slope  regimes  of  a  tw'o-component 
power-law  spectrum  of  plasma  density 
fluctuations.  The  means  we  have  used  to 
determine  the  freezing  break  is  to  observe 
the  measured  spectrum  of  plasma  density 
fluctuations  (or  of  radio  waves  passing 
through  those  fluctuations).  Our  analysis 
uses  three  different  data  sets  collected  in 
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the  naturally  disturbed  equatorial  iono¬ 
sphere.  These  are:  (1)  RF  phase  scintilla¬ 
tion  measurements.  (2)  high-resolution 
satellite  in  situ  measurements  scanning  the 
structure  horizontally,  and  (3)  in  situ 
PLUMEX  rocket  observations  scanning  the 
structure  vertically.  This  combination  of 
techniques  provides  insight  into  the  inte¬ 
grated  plasma  effects  that  impact  systems, 
the  altitude  (collisionality)  dependence  of 
the  freezing  process,  and  any  spectral 
anisotropy  associated  with  freezing. 

In  the  next  section,  we  examine  these 


spectral  forms  using  weak  scatter  phase 
scintillation  data.  In  Section  3  we  compare 
these  results  to  spectra  from  a  number  of 
orbits  of  the  AE-  E  satellite.  This  analysis 
verifies  that  the  propagation-derived  esti¬ 
mates  of  horizontally  sampled  spectral 
form  agree  with  in  situ  observations.  This 
is  followed,  in  Section  4.  by  a  detailed 
analysis  of  the  height  variations  of  spec¬ 
tral  shape  sampled  vertically  by  the  DNA 
PLUMEX  rocket,  and.  in  Section  5.  by  a 
discussion  of  the  altitude  dependence  of 
the  freezing  scale.  We  summarize  in 
Section  6. 


SECTION  2 


PHASE  SCINTILLATION  DATA:  TWO-DIMENSIONAL  SPECTRA 


The  DNA  Wideband  satellite  {Fremouw  et 
al..  1978)  operated  from  spring  1976  until 
late  summer  1979.  During  those  years, 
hundreds  of  passes  of  equatorial  phase 
scintillation  were  obtained  at  Ancon.  Peru, 
in  the  Atlantic  sector  and  Kwajalein, 
Marshall  Islands,  in  the  Pacific.  At  the  time 
of  the  first  Wideband  analyses,  there  had 
been  no  systematic  study  of  ionospheric 
irregularity  spectral  shape.  The  Wideband 
experiment  (and.  in  particular,  the  phase 
scintillation  measurements)  provided  infor¬ 
mation  for  the  first  morphological  studies 
of  equatorial  irregularity  structure.  As  a 
result,  much  of  the  application  phenome¬ 
nology  modeling  of  irregularity  spectral 
shape  has  been  based  on  these  data. 

Rino  (1979)  presents  scintillation  statistics 
of  a  large  accumulation  of  the  Wideband 
equatorial  data.  Given  a  homogeneous  and 
weakly  disturbed  medium  with  one-dimen¬ 
sional  power-law  index  p|.  phase  screen 
theory  predicts  that  the  phase  spectrum  of 
a  propagated  signal  passing  through  the 
medium  will  have  index  pi+1.  Because  that 
is  an  integrated  measure  of  structure  along 
the  raypath,  the  propagation  data  are  not 
the  unambiguous  measure  of  irregularity 
structure  that  in  situ  data  are;  in  practice, 
diffraction  and  other  propagation  effects, 
and  superposed  structured  regions,  can 
distort  a  phase  spectrum.  Nevertheless,  in 
collection,  the  average  phase  spectral  form 
can  be  measured  and  related  to  that  of  the 
irregularities.  Furthermore,  the  phase 
scintillation  data  have  direct  relevance  to 
the  prediction  of  RF  communication  and 
surveillance  system  performance,  which  is 
much  of  the  original  motivation  for  the 
experimental  work. 


At  the  time  of  the  analysis  of  the  data 
presented  by  Rino  (1979),  the  generally 
two-component  nature  of  irregularity  spec¬ 
tra  had  not  yet  been  recognized,  and  single 
power-law  fits  were  made  to  the  spectra 
over  an  approximate  spatial  wavelength 
range  of  300  to  3000  m  (temporal  fre¬ 
quency  range  0.1  to  10  Hz),  as  shown  by 
Rino  (1979).  the  average  weak-scatter  (0.2 
<  S4  <  0.5)  power-law  phase  spectral  index 
for  the  equatorial  data  was  between  2.6  and 
2.9,  corresponding  to  a  one-dimensional  in 
situ  index  between  1.6  and  1.9.  As  pointed 
out  by  Livingston  and  Dabbs  (1986).  when 
the  same  data  are  analyzed  using  more 
appropriate  spectral  analysis  techniques 
than  those  used  in  the  original  analysis 
(i.e.,  the  temporal  data  are  windowed),  the 
spectral  indices  are  somewhat  steeper:  2.9 
at  Kwajalein  and  3.1  at  Ancon. 

In  order  to  look  more  carefully  at  the  data 
in  terms  of  a  freezing  break  and  multicom¬ 
ponent  spectral  form,  a  selected  portion  of 
the  Wideband  equatorial  data  were  repro¬ 
cessed.  Some  results  of  that  analysis,  as 
applied  to  a  characterization  of  the  phase 
spectrum  under  both  weak  and  strong 
scattering  conditions,  are  discussed  by 
Livingston  and  Dabbs  (1986).  A  major 
concern  when  using  phase  spectra  to  imply 
irregularity  structure  is  that  the  particular 
manifestations  of  diffraction  on  the  phase 
spectrum  are  unknown.  By  diffraction,  we 
mean  any  departure  from  the  nominal 
linear  dependence  of  signal  phase  on 
carrier  frequency.  Unless  the  diffraction 
effects  can  be  specifically  identified,  any 
conclusions  regarding  phase  spectral  shape 
should  be  based  on  weak  scatter  data  only. 
Under  those  conditions,  phase  screen 
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theory  predicts  that  the  phase  spectrum  is 
an  integrated  mapping  of  the  in  situ 
structure,  as  noted  above.  Our  purpose  is 
to  obtain  representative  values  of  the  phase 
spectral  slopes  and  break  frequency,  from 
which  the  one-dimensional  in  situ  irregular¬ 
ity  spectral  form  can  be  implied. 

A  convenient  range  of  disturbance  level  to 
work  with  which  satisfies  the  weak  scatter 
criterion,  is  an  intensity  scintillation  level 
of  0.2  <  S4  <  0.5.  The  lower  limit  is 
necessary  to  avoid  bias  of  the  results 
caused  by  noise  contamination.  The  upper 
limit  follows  from  analysis  presented  below 
that  shows  that  above  this  level  of  S4,  phase 
errors  begin  to  rapidly  increase  and  the 
phase  spectral  shape  begins  to  be  distorted. 
There  is  an  enormous  dynamic  range  of 
equatorial  irregularity  strength;  depending 
upon  the  strength,  the  scintillation  data  that 
fall  into  this  weak  scatter  disturbance  range 
may  occur  on  any  of  the  three  Wideband 
carrier  frequencies  (137,  413,  or  1239 
MHz). 

Limiting  the  analysis  to  weak  scatter  levels 
is  not  without  its  complications.  In 
particular.  Fresnel  null  structure  is  super¬ 
imposed  on  most  of  the  weak  scatter  phase 
spectra  we  observe.  Figure  1  show's  an 
example  sequence  of  simultaneous  413- 
MHz  and  1239-MHz  intensity  and  phase 
spectra  showing  Fresnel  effects. 

(Note  that  these  1239-MHz  spectra  are  not 
included  in  our  weak  scatter  results  because 
their  S4  levels  were  less  than  0.2.)  The 
arrows  indicate  the  expected  frequency  of 
occurrence  of  the  first  Fresnel  intensity 
peak  and  phase  null,  for  the  given  pass 
geometry,  an  assumed  irregularity  anisotro¬ 
py  of  50:1  and  a  layer  height  of  350  km. 
The  nulls  are  conspicuous,  even  in  these 
spectra  that  have  been  smoothed  for 
display,  and  persist  even  beyond  the  weak 


scatter  limit.  To  avoid  Fresnel  null  biases 
w'hen  extracting  the  spectral  slopes,  two 
linear  least-square  fits  are  obtained  over 
temporal  frequency  ranges  near  the  low- 
and  high-frequency  ends  of  the  phase 
spectrum.  Taking  into  consideration  the 
average  Fresnel  null  frequencies  at  413  and 
1239  MHz  ( Livingston  and  Dabbs ,  1986), 
convenient  fit  ranges  are  0.10  to  0.56  Hz 
(-1.00  to  -0.25  log  Hz,  approximately  4500 
to  800  m  spatial  scale),  and  1.8  to  10.0  Hz 
(0.25  to  1.00  Hz,  approximately  250  to  45 
m  spatial  scale).  The  means  of  obtaining 
the  slopes  over  these  fit  ranges,  and  the 
spectral  analysis,  per  se,  are  described  by 
Livingston  and  Dabbs  (1986). 

We  chose  24  Wideband  orbits  from  Ancon 
and  19  from  Kwajalein  for  the  analysis.  All 
were  selected  because  they  show'  consistent 
scintillation  throughout  most  of  each 
pass— that  is,  at  a  variety  of  propagation 
geometries.  The  Ancon  data  are  from  the 
period  February  through  April  1977;  the 
Kwajalein  data  are  from  July  through 
September  1977.  There  is  a  distinct  shift  in 
the  season  during  which  the  irregularities 
maximize  in  the  Atlantic  and  Pacific 
sectors  (Tsunoda,  1985).  The  reprocessed 
data  are  taken  from  the  peak  scintillation 
seasons  at  each  longitude. 

The  measured  spectra  have  been  catego¬ 
rized  into  three  shape  types:  (1)  single 
power  law,  (2)  two-component  power  law 
with  a  shallow'  long  wavelength  regime 
slope  (p„)  breaking  to  a  more  steeply  sloped 
dissipation  regime  slope  (pt,),  and  (3)  the 
reverse  two-component  form.  This  catego¬ 
rization  is  somewhat  arbitrary,  and  is  not 
meant  to  imply  that  the  different  forms  are 
manifestations  of  different  irregularity 
generation  mechanisms.  Instead,  there  is  a 
continuum  of  forms  that  encompass  these 
three  types. 
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Figure  1.  Intensity  and  Phase  Spectra  for  413  and  1239  MHz  from  Ancon  Wideband  Pass.  Day  63.  1977.  Showing 
Fresnel  Effects.  The  relative  Fresnel  scale  locations  are  noted  with  arrows. 


When  the  two  measured  slopes  are  within 
20rr  of  one  another,  they  are  considered 
single-component  power-law  spectra.  As 
Table  I  shows,  these  comprise  approxi¬ 
mately  30rr  of  the  overall  spectral 
population.  Spectra  of  the  third  category 
are  those  for  which  p;i  is  significantly 
steeper  than  pi,.  Although  this  form 
accounts  for  nearly  30rc  of  the  spectra  at 
Ancon  (and  nearly  40rr  at  Kwajalein).  we 
will  not  consider  those  spectra  here.  Our 
justification  is  that  this  spectral  form  is 
almost  never  observed  in  situ  at  the  equator 
(as  will  be  shown  in  Section  3).  Most  likely, 
the  phase  spectra  of  this  type  are  a  result 
of  integration  through  more  than  just  a 
homogeneous  F  layer;  simulations  have 
shown  that  superposition  of  the  phase 
energy  from  a  high-density  E  region  can 
readily  produce  such  a  spectral  form. 

The  average  measured  Wideband  phase 
spectral  slope  data  are  summarized  in 
Table  I.  as  a  function  of  receiving  site  and 
carrier  frequency.  The  slope  distributions 
from  which  these  averages  are  derived  are 
shown  in  Figure  2  (single  power-law 
population).  Figure  3  (Ancon  two-compo¬ 
nent  population),  and  Figure  4  (Kwajalein 
two-component  population).  In  order  to 
simplify  comparison  with  the  in  situ  data, 
the  average  indices  in  Table  1  are  given  in 
terms  of  the  equivalent  one-dimensional 
form,  as  would  be  measured  by  a  satellite 
or  rocket  probe.  As  previously  mentioned, 
a  difference  of  unity  between  the  slopes  of 
phase  and  in  situ  spectra  is  predicted  by 
phase  screen  propagation  theory  at  these 
weak  scatter  levels.  The  widths  of  the 
distributions  in  Figures  2  through  4  are  also 
summarized  in  Table  1  in  terms  of  their 
I Oth  and  90th  percentile  one-dimensional 
equivalents. 

There  are  several  important  points  to  be 
made  from  Table  I,  in  addition  to  the 


measured  spectral  index  values.  The  first 
is  the  similarity  of  spectral  slopes  (both  one 
and  two  component)  at  the  different  carrier 
frequencies.  Propagation  theory  predicts 
that  although  the  energy  in  the  phase 
spectra  will  vary  from  carrier  to  carrier,  the 
spectra!  shape  will  not.  The  agreement 
between  the  slope  values  measured  at 
different  carriers  (e.g..  Ancon  137/413/ 
1239  MHz)  assures  that  the  fits  are  being 
made  without  the  influence  of  Fresnel  null 
effects.  The  agreement  between  spectral 
shapes  at  different  carriers  also  says  that 
the  spectral  form  is  not  dependent  upon 
disturbance  level.  The  irregularity  structure 
producing  the  weak  scatter  scintillation  at 
1239  MHz  is  several  times  stronger  than 
that  causing  the  same  disturbance  at  413 
MHz.  yet  their  spectral  forms  are  nearly 
identical. 

It  is  also  instructive  to  look  at  the  slope 
values  in  Table  1  relative  to  previously 
published  results.  We  have  already  men¬ 
tioned  the  Rino  (1979)  results  and  the 
limitations  of  their  applicability.  Figure  12 
of  Livingston  and  Dabbs  (1986)  summarizes 
the  average  weak  scatter  spectral  slope  of 
the  Ancon  and  Kwajalein  data,  assuming 
that  all  of  the  spectra  are  single  power  law 
in  form.  For  Ancon,  their  overall  power-law 
slope  is  close  to  that  of  the  single  power-law 
populations  in  Table  1  at  the  three  carriers. 
Furthermore,  the  average  of  the  Ancon 
two-component  indices  in  Table  1  (which 
is  a  good  approximation  to  what  would  be 
obtained  from  a  single  power-law  fit  over 
the  same  spectra)  is  also  close  to  the  same 
value.  However,  for  the  Kwajalein  data. 
Livingston  and  Dabbs  (1986)  reported  a 
significantly  lower  average  spectral  slope 
(about  1.9  in  situ)  from  single  power-law 
fits  over  all  of  the  spectra,  relative  to 
Ancon.  Until  recently,  this  has  been 
generally  considered  to  be  a  result  of 
different  (but  unidentified)  structuring 
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Table  1.  Summary  of  phase  spectral  characteristics. 


Wideband  1977-1978 


Kwajalein  (413  MHs,  254  samples) 


Power  Law  Type 

Population 

(%) 

Slope 

Slope  Range 
10th-90th  Percentile 

one 

30 

2.2 

1.5  to  2.8 

two 

33 

1.3 

0.6  to  2.0 

Pa  <  Pa 

3.0 

2.0  to  4.0 

two 

37 

2.7 

1.9  to  3.4 

Pa  >  Pa 

0.9 

0.2  to  1.6 

Break  Scale 
(m) 


Break  Range 
10th-90th  Percentile 


Power  Law  Type 

Population 

(%) 

Slope 

one 

25 

2.0 

two 

34 

1.4 

Pa  <  Pa 

3.2 

two 

41 

2.5 

Pa  >  Pa 

1.0 

Kwajalein  (1239  MHz,  76  samples) 


Slope  Range 
10th-90th  Percentile 


0.7  to  2.0 
2.3  to  4.1 


1.8  to  3.2 
0.3  to  1.7 


Break  Scale 

(m) 


Break  Range 
10th-90th  Percentile 


200  to  1700 


Ancon  (137  MHz,  110  samples) 


Power  Law  Type 

Population 

(%) 

Slope 

Slope  Range 
10th-90th  Percentile 

one 

32 

2.3 

1.8  to  2.9 

two 

38 

1.3 

0.5  to  2.2 

P«  <  Pa 

3.1 

2.1  to  4.2 

two 

30 

2.5 

1.6  to  3.3 

p«  >  Pa 

0.9 

0.2  to  1.6 

Break  Range 
10th-90th  Percentile 


Ancon  (413  MHz,  317  samples) 


Power  Law  Type 

Population 

(%) 

Slope 

Slope  Range 
10th-90th  Percentile 

one 

26 

2.1 

1.5  to  2.7 

two 

46 

1.2 

0.3  to  3.4 

Pa  <  Pa 

3.2 

2.1  to  4.2 

two 

28 

2.6 

1.8  to  3.4 

Pa  >  Pa 

1.1 

0.4  to  1.8 

Break  Scale 
(m) 


Break  Range 
10th-90th  Percentile 


Ancon  (1239  MHs,  103  samples) 


Power  Law  Type 

Population 

(%) 

Slope 

Slope  Range 
10th-90th  Percentile 

one 

19 

2.0 

1.2  to  2.8 

two 

62 

1.1 

0.3  to  2.0 

Pa  <  Pa 

3.3 

2.2  to  4.3 

two 

18 

2.7 

2.0  to  3.5 

P«  >  Pa 

0.7 

-  0.4  to  1.9 

Break  Scale 
(m) 


Break  Range 
10th-90th  Percentile 


01  23456701  234567 
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Figure  2.  Distributions  of  Single  Power-Law  Spectral  Indices  for  413  and  1239  MHz  in 
Weak  Scatter  at  Kwajalein  and  Ancon. 
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Pigure  3.  Distributions  of  the  Pa  and  Pt1  Phase  Spectral  Indices  for  413  and  1239  MHz, 
Weak  Scatter  Conditions  at  Kwajalein. 
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Figure  4.  Distributions  of  the  Pa  and  Pi,  Phase  Spectral  Indices  for  413  and  1239  MHz. 
Weak  Scatter  Conditions  at  Ancon. 
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processes  at  the  two  longitudes.  However, 
Table  1  shows  that  the  Kwajalein  single¬ 
component  and  “normal"  two-component 
slopes  are  very  similar  to  those  at  Ancon. 
We  conclude  that  the  longitude  differences 
reported  in  Livingston  and  Dabbs  (1986) 
arise  from  the  inclusion  of  the  “reverse" 
two-component  form  spectra  in  their 
analysis  population.  The  significantly  high¬ 
er  population  of  these  spectra  at  Kwajalein 
exaggerates  the  spectral  shallowing  there, 
relative  to  Ancon. 

Livingston  ct  al.  (1981).  also  made  a  study 
of  the  spectral  shape  of  a  limited  segment 
of  Kwajalein  Wideband  data  over  a  wide 
range  of  scatter  levels.  A  conclusion  of  that 
report  was  that  the  irregularity  spectral 
slope  shallowed  with  increasing  irregularity 
strength.  It  now  understood  that  phase 
spectra  shallow  during  strong  scatter 
conditions  because  of  diffraction  ( Living¬ 
ston  and  Dabbs.  1986).  We  note  that  the 
weak  scatter  portion  of  the  Livingston  et  al. 
(1981)  data,  which  were  fit  as  single-com¬ 
ponent  power  law  spectra,  exhibit  slopes 
values  close  to  those  we  observe  here. 

In  addition  to  the  slopes,  the  spatial  scale 
at  which  the  two-component  phase  spec¬ 
trum  breaks  (the  freezing  scale)  can  be 
determined.  In  the  propagation  case,  a 
break  in  the  irregularity  spectrum  appears 
in  the  phase  spectrum  at  a  temporal 
frequency  dictated  by  the  line-of-sight  scan 
velocity.  This,  in  turn,  depends  upon  the 
irregularity  altitude,  anisotropy,  and  irregu¬ 
larity  motion.  The  orbital  inclination  of  the 
Wideband  satellite  is  such  that  for  the 
premidnight  pass  data  used  here,  the 
propagation  path  scans  nearly  along  the 
magnetic  meridian.  For  F-region  irregulari¬ 
ties  that  are  highly  extended  along  the 


magnetic  field  (50:1  or  larger  anisotropy), 
the  effective  scan  velocity  caused  by  the 
satellite  motion  depends  upon  propagation 
geometry,  but  has  an  average  value  of  ~ 
0.45  km/s.  At  the  local  time  and  latitude 
of  the  observations.  F-region  polarization 
and  tidal  electric  fields  will  typically  drive 
the  irregularities  eastward  at  a  velocity  of 
~  0.25  km/s.  Thus  we  can  estimate  the 
total  effective  scan  velocity  to  be  in  the  0.5- 
to  0.7-km/s  range. 

For  each  two-component  spectrum  in  the 
population  (both  normal  and  reverse 
forms),  the  break  frequency  is  determined 
as  the  intersection  of  the  p;,  and  pi,  linear 
fit  segments.  The  distribution  of  the  break 
frequencies  is  shown  in  Figure  5  for  413 
and  1239  MHz  at  Ancon  ami  Kwajalein. 
The  distributions  are  broad,  which  is 
primarily  a  result  of  the  range  of  sampled 
propagation  geometries  and  effective  scan 
velocities.  If  there  were  a  means  of 
accurately  apply  ing  corrections  for  the  scan 
velocity  so  that  the  data  could  be  sorted  in 
spatial  scale,  the  distributions  would  be 
much  narrower. 

Despite  the  broadening,  the  distributions  in 
Figure  5  all  peak  near  the  same  frequency. 
Because  the  temporal  frequency  at  which 
the  translated  break  occurs  is  independent 
of  carrier  frequency,  this  lends  credence  to 
the  interpretation  that  the  origin  of  the 
break  is  in  situ.  Choosing  an  overall 
average  break  frequency  of  1.1  Hz.  and  a 
representative  scan  velocity  of  0.60  km/s. 
our  estimate  of  the  irregularity  break  scale 
is  near  0.55  km.  This  is  in  good  agreement 
with  the  irregularity  spectrum  break  scale 
observed  by  satellites  and  rockets,  as  will 
be  reviewed  in  the  following  sections. 


PERCENTAGE 


Figure  5.  Distributions  of  Break  Frequency  for  413  and  1239  MHz  Weak  Scatter 
Measured  at  Kwajalein  and  Ancon. 
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SECTION  3 


AE-E  DATA:  HORIZONTAL,  ONE-DIMENSIONAL  SCANS 


The  Atmospheric  Explorer-E  (AE-E)  satel¬ 
lite  mission  operated  from  late  1976 
through  micll 979,  providing  in  situ  mea¬ 
surements  of  the  equatorial  ionosphere. 
Among  the  instruments  that  the  spacecraft 
carried  was  a  retarding  potential  analyzer 
(RPA)  from  which  high-resolution  mea¬ 
surements  of  plasma  density  structure  were 
obtained.  A  detailed  description  of  the  RPA 
and  its  associated  instrumentation  can  be 
found  in  work  by  Hanson  and  Heelis  (1975). 

Of  interest  to  us  here  are  two  periods  of 
time  when  AE-E  was  operating  at  altitudes 
well  within  the  equatorial  F  layer.  These 
two  periods  were  in  autumn  1978  and 
w  inter  1979-1980,  when  the  satellite  was  at 
altitudes  of  375  and  435  km.  respectively. 
This  first  period  was  well  onto  the  upswing 
of  the  solar  cycle:  the  latter  period  w  as  near 
sunspot  maximum.  During  each  period, 
many  orbits  of  data  were  obtained.  Of 
these,  there  were  approximately  30  orbits 
when  the  RPA  instrument  was  on  and  the 
ionosphere  that  it  sampled  was  significantly 
structured  to  merit  detailed  analysis. 
Because  at  low  latitudes  the  F  region 
strictures  only  during  the  local  evening  and 
night  hours  when  F-layer  bubbles  can 
occur,  the  disturbed  orbits  all  correspond 
to  local  evening.  The  previously  mentioned 
seasonal  dependence  of  irregularity  occur¬ 
rence  effectively  limits  the  selected  1978 
data  to  Pacific  longitudes,  and  the 
1979-1980  data  to  Atlantic  longitudes. 

The  RPA  instrument  aboard  AE-E  did  not 
operate  continuously,  but  rather  collected 
data  for  three  out  of  each  six  seconds. 
Figure  6  shows  a  series  of  these  samples 
with  respect  to  the  continuous  low- 


resolution  ion  density  measured  by  the  drift 
meter.  For  the  orbits  of  interest,  the  RPA 
data  spans  that  exhibited  structure  (and 
some  degree  of  stationarity  in  that  struc¬ 
ture)  were  selected  for  processing.  The  only 
structures  that  w'ere  specifically  avoided 
were  those  obviously  dominated  by  bottom- 
side  spread  (BSS)  irregularities;  these  were 
found  primarily  in  the  375-km  data.  For 
each  RPA  sample,  background  parameters 
(e.g.,  average  ion  density,  density  gradient) 
have  been  extracted,  and  the  power 
spectral  density  of  the  plasma  structure 
computed. 

As  an  overview  of  the  background  condi¬ 
tions,  Figure  7a  shows  the  variation  of 
mean  ion  density  with  dip  latitude:  each 
point  corresponds  to  an  RPA  record  for 
which  spectra  have  been  computed.  The 
figure  illustrates  that  the  background 
densities  for  the  1978  (375  -  km)  and 
1979-1980  (435-km)  measurements  are 
very  similar.  Both  show  the  same  increase 
in  ion  density  away  from  the  magnetic 
equator,  which  is  probably  a  signature  of 
the  equatorial  anomaly.  The  similarity  of 
the  mean  densities  is  not  surprising,  be¬ 
cause,  relative  to  the  F  peak,  the  obser¬ 
vation  altitudes  are  nearly  the  same, 
lonosonde  measurements  collected  at  Kwa- 
jalein  during  mid-1979  ( Tsunoda  and  Rino , 
1982)  show  a  typical  F-layer  bottomside 
altitude  of  350  km  at  spread-F  onset,  with 
a  decline  to  300  km  by  midnight.  Similar 
observations  made  in  mid- 1981.  late  in  the 
peak  of  the  solar  cycle,  show  the  bottom- 
side  altitudes  to  be  shifted  upward  by  an 
average  of  50  km.  We  would  estimate  that 
in  both  our  1978  and  1979-1980  data  sets. 
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AE-E  Satellite  in  situ  Plasma  Density  Observations.  Intervals  of  high  resolution  RPA  data  are  indicated  in  the 
top  panel  and  detailed  below. 


the  AE-E  was  very  near  the  average  peak 
of  the  F  layer  at  local  middle  evening  hours. 

The  integrated  energy  in  the  overall, 
intermediate  scale  (200-  to  2000-m)  ion 
density  structure  is  shown  in  Figure  7b.  The 
peak  deviations  in  the  anomaly  region  are 
similar,  but  at  the  dip  equator,  the  average 
energy  is  stronger  at  375  km.  Figure  7c 
illustrates  that  both  data  sets  show  similar 
and  systematic  weakening  of  the  structure 
with  local  time  as  the  irregularities, 
generated  by  plasma  bubble  passage, 
slowly  decay. 

For  our  morphological  purposes  here,  the 
data  in  Figure  7  show  that,  in  collection, 
there  are  no  dramatic  differences  between 
the  two  time/altitude  periods,  in  terms  of 
background  conditions  and  level  of  struc¬ 
turing.  This  simplifies  the  comparison  of 
the  spectral  forms  obtained  from  the  two 
data  sets;  with  no  gross  differences  between 
the  ionospheric  conditions  (background 
density  or  large-scale  gradients),  the  pri¬ 
mary  structuring  processes  should  be 
similar. 

Figure  8  shows  typical  power  spectra  of  the 
ion  density  irregularities  (these  are  from 
one  orbit  in  the  1978  series).  The  spectra 
are  computed  in  a  conventional  fashion, 
using  the  procedures  developed  for  the 
HILAT  satellite  program.  The  three-second 
RPA  records  are  edge  matched  and  cosine 
windowed,  before  computation  of  a  600- 
point  FFT.  This  provides  spectral  informa¬ 
tion  over  a  frequency  range  of  0.5  to  100.0 
Hz.  For  the  satellite  scan  velocity  of  0.77 
krn/s.  this  corresponds  to  measurement  of 
structure  energy  at  spatial  wavelengths 
between  about  80  and  4000  m.  The 
computed  spectra  are  binomially  smoothed 
and  log  decimated  to  50  points,  producing 
the  form  shown  in  Figure  8. 


One  of  the  difficulties  encountered  when 
deriving  spectral  shape  information  (slopes 
and  break  scales)  from  multicomponent 
power-law  spectra  is  that  all  practical 
automated  fitting  algorithms  appear  to  have 
biases.  When  the  low-  and  high-frequency 
portions  of  a  multicomponent  spectrum 
have  very  different  slopes,  and  the  break 
itself  is  the  dominant  spectral  feature, 
accurate  fits  are  derived.  However,  when 
the  slopes  are  similar,  or  the  spectrum 
shows  some  isolated  secondary  structure, 
fits  are  often  obtained  that,  by  eye,  do  not 
appear  to  be  optimum.  For  these  AE-E 
data,  it  is  practical  to  fit  each  spectrum 
interactively.  Breaks  in  the  spectra  are 
identified  by  eye  in  a  computer  display,  and 
least-square  slope  fits  to  the  power-law 
slope  segments  between  these  breaks  are 
determined  by  the  computer.  Figure  8 
shows  that  the  power-law  segments  deter¬ 
mined  this  way  accurately  represent  the 
measured  continuum. 

The  collected  spectral  form  data  have  been 
sorted  into  various  spectral  types.  The 
populations  for  each  type  are  determined 
using  measured  spectral  breaks  that  are 
within  the  freezing  break  “window"  (150- 
to  2500-m  spatial  scale);  as  with  the  phase 
spectra  in  Section  2,  adjoining  slope 
segments  that  are  within  20%  of  one 
another  are  considered  to  be  a  single 
power-law  segment  with  no  intervening 
break.  Table  2  shows  the  resulting 
population  of  one-,  two-,  and  three-compo¬ 
nent  spectral  types.  For  both  time  periods, 
the  majority  of  the  spectra  are  two 
component  in  shape,  with  the  expected 
form:  a  shallow  intermediate  scale  regime 
breaking  to  a  more  steeply  sloped  dissipa¬ 
tion  regime.  The  reverse  form  (i.e..  more 
shallowly  sloped  at  short  scales)  is  virtually 
absent,  unlike  the  phase  spectral  data. 
Significant  portions,  roughly  one  quarter. 


15 


AE-E 

1978 


UJ 

O 

o 


< 


uj/3*(»  6o|  —  A1ISN3Q  NOI  NV3W 


16 


Figure  7.  Background  Conditions  Measured  by  the  AE-E  Satellite  in  Terms  of  Mean  Plasma  Density  and  its  Fluctuations 
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Table  2.  Summary  of  in  situ  irregularity  spectral  characteristics. 


AE-E  1978 


375  km  Altitude,  Pacific  Sector,  272  samples 

Population 

Slope  Range 

Break  Scale 

Break  Range 

Power  Law  Type 

(%) 

Slope 

10th-90th  Percentile 

(m) 

10th-90th  Percentile 

one 

23 

1.6 

1.1  to  2.1 

— 

two 

56 

0.8 

0.1  to  1.5 

800 

500  to  1300 

p«  <  Pfc 

2.4 

1.9  to  2.9 

two 

1 

— 

Pa  >  P6 

— 

— 

three 

20 

0.9 

-  0.3  to  2.1 

■ 

1200  to  2500 

2.0 

0.1  to  1.6 

1 

240  to  550 

2.4 

1.4  to  3.4 

AE-E  1979-1980 


435  km  Altitude,  Atlantic  Sector,  146  samples 

Population 

Slope  Range 

Break  Scale 

Break  Range 

Power  Law  Type 

(%) 

Slope 

10th-90th  Percentile 

M 

10th-90th  Percentile 

one 

29 

1.9 

1.3  to  2.5 

— 

67 

0.9 

0.1  to  1.6 

720 

480  to  1100 

■SH 

3.0 

2.5  to  3.5 

two 

1 

— 

Pa  >  Pft 

— 

— 

three 

4 

— 

— 

— 

of  both  the  1978  and  1979-1980  popula¬ 
tions  are  best  fit  by  a  single  power-law 
segment.  The  biggest  difference  between 
the  distributions  measured  at  the  two 
altitudes  is  the  relatively  high  percentage 
of  three-component  spectra  in  the  1979- 
1980  data. 

Figure  9  shows  the  distributions  of  the 
measured  spectral  slopes  for  the  single 
power-law  and  two-component  AE-E  spec- 
fa!  forms.  These  are  summarized  in  Table 
2  in  terms  of  the  slope  mean  and  10th/90th 
percentile  range.  Although  most  of  the 
discussion  that  follows  is  concerned  with 


the  dominant,  two-component  spectra,  it  is 
worthwhile  to  consider  the  overall  spectral 
character  that  is  observed.  As  with  the 
phase  spectra  populations,  we  would  argue 
that  although  the  spectra  are  sorted  into 
three  population  types  in  Table  2  for 
convenience,  they  are  actually  samples  of 
the  same  population  continuum.  Instead  of 
being  the  product  of  separate  structuring 
mechanisms,  they  are  manifestations  of  the 
same  process  encountered  under  different 
conditions.  For  both  altitude  periods,  the 
average  of  the  two-component  spectral 
indices  is  nearly  exactly  the  same  as  the 
measured  single  power-law  index:  this  is 
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because  the  average  of  the  indices  is  a  good 
approximation  to  what  would  be  obtained 
using  a  single  power-law  fit  over  the  entire 
spectrum.  The  same  is  true  for  the  375-km 
altitude,  three-component  indices. 

There  are  few  data  with  which  to  compare 
the  slope  information  in  Table  2,  other  than 
the  phase  data  presented  in  Section  2. 
Examples  of  in  situ  irregularity  spectra 
measured  by  AE-E  are  presented  in  a  few 
papers  (e.g.,  McClure et  al.,  1977)  but  there 
have  been  no  systematic  studies  of 
ionospheric  irregularity  spectra  in  the  200- 
to  2500-m  spatial  wavelength  “freezing 
window.”  Livingston  et  al.  (1981)  present 
spectral  results  from  low-resolution  AE-E 
data,  based  upon  single  power-law  fits 
made  at  longer  wavelengths  (750  to 
7500-m).  The  physics  controlling  the 
spectral  shape  at  those  longer  wavelengths 
appears  to  be  quite  different  from  that 
controlling  the  intermediate  scale  structure, 
and  comparisons  between  these  data  and 
those  of  Livingston  et  al.  (1981)  are  not 
appropriate. 

The  statistical  summary  of  the  irregularity 
spectral  slopes  in  Table  2  shows  that 
despite  the  similarity  of  background 
conditions  and  structuring  levels,  the  two 
different  year/altitude/-  longitude  sets  show 
some  differences.  The  435-km  single 
power-law  slopes  are  considerably  steeper 
than  those  at  375  km.  The  two-component 
slope  values  suggest  that  this  difference 
arises  from  a  steeper  spectrum  at  short 
spatial  wavelengths;  there  is  very  little 
difference  between  the  slopes  at  wave¬ 
lengths  larger  than  the  freezing  scale. 

It  is  difficult  to  determine  any  specific 
reason  for  these  differences  because  of  the 
separation  of  the  data  sets  in  solar  cycle, 
longitude  (through  magnetic  field  configu¬ 
ration),  and  altitude.  Any  of  these  factors 


can  affect  the  growth  and  dissipation 
mechanisms  that  determine  the  details  of 
the  short-wavelength  portion  of  the  spectra. 


The  measured  distributions  of  the  two- 
component  spectral  break  wavelength  are 
shown  in  Figure  10,  and  are  summarized 
in  Table  2.  For  each  spectrum  in  the 
two-component  populations,  the  break 
wavelength  data  have  been  converted  to  an 
equivalent  cross-field  dimension,  based 
upon  the  satellite  orbit  and  the  DGRF-75/80 
magnetic  field  model.  It  is  interesting  to 
note  from  Figure  10  that  the  distribution  of 
break  scales  is  relatively  uniform  over  the 
standard  deviation  range  shown  in  Table  2. 

Unlike  the  spectral  slopes,  which  differ 
between  the  two  different  altitude  sets, 
there  is  only  a  small  difference  between  the 
average  freezing  break  wavelengths  at  375 
and  435  km.  The  small  change  that  is  seen 
is  attributable  to  the  occurrence  of  more 
larger-scale  (  >  1100  m)  breaks  in  the 
375-km  data,  as  the  distributions  in  Figure 
10  illustrate. 

Although  the  specific  mechanisms  for 
control  of  the  freezing  break  are  not 
completely  understood,  it  is  likely  that 
collisionality  plays  a  significant  role.  The 
neutral  atmospheric  density,  and  therefore 
the  collision  frequency,  decreases  exponen¬ 
tially  with  altitude.  In  the  collision- 
dominated  F  region,  the  irregularity  growth 
rate  is  inversely  proportional  to  collision 
frequency;  one  would  expect,  to  first  order, 
that  this  would  result  in  an  altitude 
dependence  of  the  freezing  scale  wave¬ 
length.  It  is  therefore  somewhat  surprising 
to  find  such  close  agreement  in  break  scales 
measured  at  375  and  435  km.  Accordingly, 
we  have  considered  the  collision  conditions 
for  the  two  data  sets  in  more  detail. 
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The  average  ion-neutral  collision  frequency 
for  the  two  time  periods  of  AE-E  data  has 
been  estimated  from  measured  FI 0.7  solar 
flux.  Kelly  ct  ul.  (1983)  present  the  daily 
mean  exospheric  temperature  versus  solar 
flux  for  a  solar  cycle  of  Chatanika  radar 
observations.  A  fit  to  their  data  gives  the 
expression 

f*  =  685.  +  2.54  *F,„.7 

For  the  AE-E  data  periods  in  1978  and 
1979-1980.  the  average  estimated  exo¬ 
spheric  temperatures  are  1041  K  for  the 
Pacific  sector  data  and  1180  K  for  the 
Atlantic  sector  data.  Figure  1 1  shows  the 
altitude  variations  of  the  ion-neutral 
collision  frequencies  based  on  the  Jaccia 
neutral  density  model  and  using  the 
momentum  transfer  cross-section  coeffi¬ 
cients  after  Scluink  and  Walker  (1973).  As 
the  figure  shows,  the  intervening  year  of 
increasing  solar  flux  reduces  the  collision 
frequency  change  between  375  and  435  km. 
However,  even  the  change  from  0.20  to 
0.12  Hz  should  produce  a  measurable 
increase  in  the  freezing  break  between  the 
lower-  and  higher-altitude  data  periods;  the 
data  do  not  show  such  a  change. 

There  are  a  number  of  reasons  that  this 
first-order  application  of  collision  control 
of  the  freezing  break  might  not  apply  to  the 


AE-E  data.  The  control  of  the  freezing 
break  may  be  much  more  complex  than  we 
have  suggested,  but  there  also  may  be  a 
relatively  simple  explanation.  The  basic 
assumption  used  above  is  that  the  freezing 
break  is  set  by  a  single,  constant  driver 
(gravity)  under  conditions  in  the  irregular¬ 
ity  locale.  The  PLUMEX  rocket  data  that 
will  be  presented  in  Section  4  support  this 
assumption,  except  that  the  freezing  break 
in  those  data  appear  to  be  set  by  local 
conditions  at  the  time  and  place  that  the 
irregularities  were  created.  As  is  well 
known,  the  dynamics  of  growth-phase 
equatorial  bubbles  can  rapidly  transport 
kilometer-scale  structure  over  a  large 
altitude  range.  Furthermore,  any  dominant 
kilometer-scale  structure  should  map  along 
field  lines.  Thus,  AE-E  in  an  orbit  that 
traverses  approximately  ±20  degrees  of 
dip  latitude,  will  sample  a  complicated 
composite  of  freezing  breaks  that  reflect  a 
wide  variety  of  irregularity  generation 
conditions. 

In  contrast  to  this  broad  time  and  space 
coverage  of  AE-E.  in  the  next  section  we 
review  data  from  the  PLUMEX  II  rocket. 
These  data  provide  a  means  to  obtain 
spectral  shape  under  known  local  condi¬ 
tions  and  over  a  short  enough  period  in  time 
that  any  drivers  are  constant. 
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COLLISION  FREQUENCY  —  Hz 


Figure  11.  Collision  Frequency  Profiles  Estimated  on  the  Basis  of  the  Measured  F10.7 
cm  Flux  and  the  Jacchia  Neutral  Atmosphere  Model. 


SECTION  4 


PLUMEX  DATA;  VERTICAL,  ONE-DIMENSIONAL  SCANS 


The  DNA  PLUMEX  II  rocket  was  launched 
in  July  1979  from  Kwajalein  Atoll  into  fully 
developed,  early  decay  phase  ionospheric 
F-region  structure.  Its  trajectory  followed 
the  central  stem  of  an  equatorial  plume  that 
extended  from  the  bottomside  through  the 
topside  of  the  F  layer.  During  its  upleg 
flight,  the  rocket  passed  through  a  highly 
structured  F  region  and  then  a  distinct 
topside  depletion.  On  the  downleg  flight, 
the  rocket  again  penetrated  the  depletion, 
at  a  slightly  more  northerly  latitude,  and 
encountered  F-region  structure  similar  to 
that  seen  on  the  upleg. 

The  PLUMEX  data  have  been  discussed  in 
several  papers  and  reports.  Most  notably. 
Kcllvv  et  ul.  (1982)  studied  short-scale 
structure,  and  Rina  et  ul.  (1981 )  considered 
intermediate-scale  (100-  to  2000-rn)  irregu¬ 
larities.  Rino  et  ul.  (1981)  performed 
spectral  analysis  on  a  broad-altitude 
segment  of  the  F-layer  electron  density 
profile,  and  found  a  two-component 
power-law  form,  with  a  break  near  800  m. 
This  is  very  similar  to  the  statistical 
characterization  obtained  from  the  AE-E 
data  presented  in  Section  3.  The  implica¬ 
tion  of  this  agreement  is  that  vertical  and 
horizontal  in  situ  scans  (both  cross  field) 
are  essentially  equivalent,  (i.e.,  there  is 
little  anisotropy  in  the  plane  perpendicular 
to  the  magnetic  field)  and  that  the 
PLUMEX  data  provide  a  representative 
view  of  the  altitude  dependencies  of 
equatorial  irregularity  processes  at  interme¬ 
diate  scale  sizes. 


dominant  irregularity  structure  in  the  F 
region.  With  the  emergence  of  theories  to 
explain  multicomponent  spectra,  we  have 
returned  to  the  PLUMEX  data  to  make  a 
closer  examination  of  the  irregularity 
spectral  forms  that  were  observed.  The 
PLUMEX  experiment  was  unique  in  that  it 
provided  a  means  to  examine  altitude 
dependencies  in  the  irregularity  processes, 
a  dimension  not  generally  available  from 
satellites.  In  particular,  it  is  possible  to 
determine  the  altitude  dependence  of  the 
break  point. 

We  first  consider  the  overall  geometry  of 
the  PLUMEX  experiment,  which  has  an 
important  bearing  on  the  upleg/downleg 
spectral  comparisons.  Figure  12  shows  the 
rocket  trajectory  (solid  curve)  and  magnetic 
field  lines  (DGRF75)  at  three  altitudes, 
traced  from  the  downleg  to  the  upleg  and 
on  southward  to  the  magnetic  equator.  The 
altitudes  of  the  field  lines  at  their 
intersections  with  the  trajectory  and  the 
magnetic  equator  are  noted.  In  actuality, 
the  rocket  trajectory  was  not  in  the 
magnetic  meridian  (8.5  degrees  azimuth) 
as  Figure  12  implies;  instead  it  was 
launched  at  about  25  degrees  azimuth, 
which  allowed  interception  of  the  eastward 
drifting  plume  during  both  the  upleg  and 
downleg  flights.  Several  features  of  the 
macroscale  structure  encountered  during 
the  upleg  were  also  observed  during  the 
downleg. 

Figure  13  shows  the  plasma  density  profiles 
obtained  during  the  upleg  and  downleg. 
The  topside  bubble  (at  475  km  on  the 
upleg)  was  encountered  at  about  440  km  on 
the  downleg;  the  same  structured  sublayer 


Because  it  was  computed  over  such  a  large 
altitude  range,  the  Rino  et  ul.  (1981) 
spectral  result  provided  only  a  look  at  the 
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Figure  12.  Geometry  of  the  PLUMEX  I  Experiment.  Showing  the  Rocket  Trajectory 
and  Magnetic  Field.  Altitudes  of  the  field  line  intersections  with  the 
trajectory  and  the  magnetic  equator  are  as  indicated. 
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Figure  13.  Plasma  Density  Profile  Measured  by  the  PLUMEX  Rocket  on  the  Upleg  and 
Downleg  Respectively. 
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seen  between  355  and  410  km  upleg.  is  seen 
between  320  and  at  375  km  on  the  downleg. 
The  altitude  differences  between  the 
upleg/downleg  macroscale  plasma  features 
coincide  well  with  the  altitude  variation  of 
the  magnetic  field  lines  shown  in  Figure  12. 
The  implication  is  that  the  macroscale 
plasma  features  in  the  middle  and  upper  F 
layer  map  from  upleg  to  downleg  along  the 
magnetic  field.  The  profiles  also  suggest 
that  local  aeronomy,  rather  than  mapping, 
controls  the  layer  shape  at  low  altitudes. 
Using  AI.TAIR  incoherent-scatter  radar 
measurements  made  during  a  scan  in  the 
magnetic  meridian.  Tsitnoda  (1980)  showed 
similar  field  alignment  of  macroscale 
plasma  structure  at  high  altitudes,  and  the 
control  of  the  bottomside  F  layer  by 
recombination  chemistry.  In  our  case,  this 
handoff  between  mapping  and  aeronomy 
processes  occurs  at  an  altitude  of  about  275 
km.  Similar  conclusions  were  drawn  by 
Kelley  et  al.  (1986)  in  the  CONDOR  experi¬ 
ments. 

The  macroscale  similarities  between  the 
upleg  and  downleg  can  be  seen  more 
conveniently  in  Figure  14.  where  the  mean 
density  and  gradient  scale  are  plotted 
versus  altitude.  Both  parameters  are  com¬ 
puted  over  a  10-km  interval  in  altitude.  In 
the  figure,  the  downleg  data  have  been 
shifted  upward  in  altitude  by  the  amount 
determined  by  the  field  line  trace,  referenc¬ 
ing  them  to  the  same  magnetic  altitude  as 
the  upleg. 

Point  by  point,  there  are  differences 
between  the  profiles,  but  the  large-scale 
plasma  sublayers  and  the  gradient  scales 
are  quite  similar,  overall. 

On  the  basis  that  the  bulk  plasma  features 
map  from  the  upleg  to  the  downleg,  one 
might  expect  that  the  intermediate-scale 
structure  would  also  map.  This  does  not 


seem  to  be  the  case.  Figure  15  shows  that 
irregularity  structure  on  the  upleg  is 
concentrated  on  the  positive  gradients  of 
the  profile;  this  would  imply  that  the 
positive  gradients  were  unstable  to  convec¬ 
tive  structuring  prior  to,  or  during,  the 
upleg  flight.  In  contrast,  on  the  downleg 
profile  there  is  no  similar  correlation 
between  gradient  sign  and  irregularity 
strength.  Figure  15  specifically  compares 
the  altitude  distribution  of  the  intermediate- 
scale  structure,  after  mapping  the  downleg 
statistics  to  the  upleg  altitude.  The  first 
parameter  shown  is  the  logarithm  of  the 
standard  dev  iation  in  electron  density.  oN. 
The  second  is  the  logarithm  of  the 
normalized  deviation.  oN/N.  Because  we 
are  concerned  specifically  with  the  struc¬ 
ture  in  the  freezing  scale  range,  the 
parameters  are  obtained  by  integrating  the 
spectral  energy  over  a  2000-  to  200-m  wave¬ 
length  window. 

Because  of  the  power-law  nature  of  the 
irregularity  spectrum,  both  parameters  in 
Figure  15  are  dominated  by  irregularities 
at  the  largest  scales  in  this  range,  and  these 
should  effectively  map  along  the  magnetic 
field.  However  Figure  15  shows  that  the 
structure  maxima  do  not  correspond, 
indicating  that  the  1-  to  2-  km  wavelength 
density  irregularities  do  not  map  along  the 
field  from  upleg  to  downleg.  The  most 
likely  explanation  is  that  the  upleg  and 
downleg  trajectories  penetrated  slightly 
different  parts  of  the  plume— the  upleg 
sampling  the  central  stem,  and  the  downleg 
sampling  an  edge  where  horizontal  gradi¬ 
ents  could  be  important.  An  alternative,  but 
probably  less  likely  explanation  in  this 
case,  is  that  irregularities  with  a  real 
frequency  component  can  transport  them¬ 
selves  from  the  gradients  on  which  they 
were  produced  (Sperling  and  Glassmun. 
1985).  For  our  purposes  here,  however,  the 
precise  altitude  distribution  of  irregularities 
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14.  Macroscale  Plasma  Density  Features  Measured  bv  the  PLUMEX  Rocket 


is  less  important  than  the  fact  that  the 
perturbation  levels  of  the  upleg  and 
downleg  data  are  similar:  this  makes  the 
direct  comparison  of  their  higher-order 
statistics  more  straightforward. 

In  this  investigation,  we  have  attempted  to 
maximize  the  altitude  resolution  of  the 
spectral  analysis,  in  order  to  best  observe 
height  dependencies  of  the  break  scales. 
This  is  easy  to  implement,  but  the  results 
can  be  difficult  to  interpret  in  terms  of 
trends  in  spectral  shape.  The  best  resolu¬ 
tion  that  is  practical  is  about  10  km  in 
altitude,  which  is  the  minimum  span 
suitable  for  spectral  analysis.  At  this 
resolution,  however,  there  is  significant 
variability  of  spectral  form  from  sample  to 
sample,  even  when  the  structuring  process 
is  clearly  dependent  on  altitude.  For  the 
analysis  here,  we  have  tried  to  simplify  the 
spectral  shape  interpretation  by  selecting  a 
series  of  20-km  spans  of  data  that:  (1) 
appear  to  be  stationary,  and  (2)  are  regions 
over  which  the  irregularity  generation 
process  is  likely  to  be  consistent.  For 
example  the  spans  selected  for  the  upleg 
are  the  positive  and  negative  gradients  of 
the  sublayers,  as  noted  in  Figure  14.  This 
selection  of  specific  spans  assures  that  if 
the  break  scale  depends  on  either  the  driver 
or  the  disturbance  level,  the  difference  in 
break  scale  between  those  spans  will  not  be 
obscured.  Within  each  20-km  span,  spectra 
are  computed  from  10-km.  overlapped, 
altitude  ranges  (the  top,  bottom,  and 
middle  10  km  of  the  data).  By  using  three 
individual  spectra,  we  can  convey  some 
qualitative  measure  of  the  spectral  variabil¬ 
ity  (or  consistency)  which  can  occur  over 
each  narrow  (20-km)  altitude  range. 

Figure  16  shows  data  from  one  low-altitude 
portion  of  the  upleg.  Because  spectral 
shape  is  our  primary  concern,  the  three 
individual  (10-km  overlapped)  spectra  are 


averaged  in  log  power  to  obtain  the  mean 
shape  (for  the  overall  20-km  span)  shown 
as  a  dashed  line.  Note  that  the  mean  shape 
is  not  the  same  as  the  mean  spectral 
density,  which  would  follow  the  top  of  the 
envelope  formed  by  the  three  curves.  A  fit 
is  applied  to  the  mean  spectral  shape  in 
Figure  16  using  the  same  interactive 
procedure  as  was  used  with  the  AE-E  data. 
That  is.  break  points  are  identified  by  eye. 
and  slopes  fits  are  made  by  the  computer 
in  a  least-squares  sense. 

The  spectra  in  Figure  16  are  from  the 
positiv  e  gradient  region  on  the  upleg  F  layer 
bottomside  (260-  to  310-km).  All  three 
altitude  intervals  are  consistent  in  shape, 
both  individually  and  on  average,  with  a 
spectral  break  near  300-m  wavelength. 

Between  the  bottom  and  the  top  of  the 
20-km  altitude  span,  there  is  a  large 
difference  in  the  fluctuation  energy  at  long 
wavelengths,  but  at  scales  shorter  than  the 
break  the  slopes  and  energies  are  almost 
identical. 

There  are  six  other  stationary  data  spans 
for  the  upleg  flight  (see  Figure  13)  in 
addition  to  that  shown  in  Figure  16.  The 
density  data  and  spectra  for  four  of  these 
are  shown  in  Figure  17.  Both  positive  and 
negative  gradient  regions  are  represented 
among  the  segments;  the  clear  differenti¬ 
ation  between  the  amount  of  density 
structure  suggests  that  these  have  been 
recently  driven  and  nondriven,  respective¬ 
ly.  Also  clear  from  Figure  17  is  that  there 
is  a  variation  of  the  freezing  break  with 
altitude.  Of  the  two  other  data  spans  not 
shown  in  the  figure,  one  (segment  b  in 
Figure  13)  is  very  similar  in  form  to  Figure 
15.  The  other  (segment  g  in  Figure  13)  is 
on  the  topside.  It,  too.  shows  a  spectral 
break,  but  the  spectral  slopes  are  similar 
enough  that  the  overall  form  can  also  be 
interpreted  as  single  power  law. 
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Figure  16.  A  20  km  Altitude  Segment  of  the  PLUMEX  I  Upleg  Electron  Density  Profile 
and  its  Spectral  Form.  Shown  are  the  spectra  from  three  overlapped  10-km 
spectra  (dotted),  the  average  spectral  shape  (dashed),  and  the  least  squares 
fit  to  the  average  shape. 
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Six  segments  of  data  from  the  downleg  also 
have  been  analyzed.  Telemetry  data  gaps 
limit  these  useful  spans  to  altitudes  below 
380  km.  as  shown  in  Figure  13.  As  noted 
above,  the  macroscale  features,  but  not  the 
kilometer-scale  structure,  appear  to  map 
from  the  upleg  to  the  downleg  along  the 
magnetic  field.  It  can  be  safely  assumed 
that  the  irregularity  structures  seen  upleg 
and  downleg  were  generated  at  approxi¬ 
mately  the  same  time,  and  through  a 
similar  process.  The  downleg  structure  is 
only  slightly  weaker  in  strength  and  is 
sampled  about  ten  minutes  after  the  upleg. 
so  there  should  be  little  diffusive  loss  in 
spectral  energy.  In  fact,  this  appears  to  be 
the  case:  the  downleg  spectra  are  similar 
to  those  in  Figure  17  in  terms  of  break 
wavelength  variation  and  slope.  The  pri¬ 
mary  difference  is  that  there  is  no  ordering 
of  the  irregularity  strength  with  gradient 
sign,  as  there  is  on  the  upleg. 

The  break  wavelength  data  from  all  the 
upleg  and  downleg  spectra  are  combined 
in  Figure  18.  The  direction  of  the  arrows 
at  each  point  indicates  which  are  from  the 
upleg  and  which  are  from  the  downleg.  In 
order  to  combine  the  data  on  a  common 
vertical  scale,  we  have  used  magnetic 
altitude— that  is.  the  height  at  which  the 
magnetic  field  line  for  each  sample 
segment  crosses  the  magnetic  equator. 


In  effect,  this  simply  means  that  the 
downleg  samples  are  shifted  upward  about 
30  km  relative  to  the  upleg  samples,  ac¬ 
cording  to  the  magnetic  field  configuration 
in  Figure  12.  The  shift  is  not  essential  to 
our  conclusions,  but  brings  the  upleg  and 
downleg  samples  into  very  close  alignment. 
This  alignment  is  the  most  obvious  feature 
of  Figure  18.  There  is  a  systematic  ordering 
of  all  but  the  topside  and  bottomside  break 
scales.  They  range  from  about  300  m  at  270 
km  altitude,  to  larger  than  a  kilometer  at 
400  km  altitude. 

There  is  also  a  change  of  spectral  slope  with 
altitude,  as  the  spectra  in  Figures  16  and 
17  imply.  Overall,  the  indices  of  the  two 
power-law  components  of  the  spectra  both 
become  more  shallow  with  increasing  al¬ 
titude.  The  average  slopes  for  the  two- 
component  spectra  whose  breaks  lie  along 
the  main  sequence  in  Figure  18  are  p;i  =  0.8 
and  pi,  =  3.3.  The  shallower  slope  is  in  good 
agreement  with  the  horizontal  scan  AE-E 
value  (Table  2):  the  dissipation  regime 
slope  is  significantly  steeper  than  the  AE-E 
average.  As  mentioned  above,  the  highest 
altitude  spectrum  is  better  characterized  as 
single  power  law  in  form.  This  is  also  true 
for  the  next  spectrum  down  in  altitude, 
where  the  break  wavelength  approaches  the 
large  scale  limit  of  the  spectrum.  For  these 
topside  samples,  the  spectral  form  ap¬ 
proaches  k"2. 
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Figure  18.  Altitude  profile  of  bread  scales  measured  by  PLUMEX 


3' 


SECTION  5 


THE  ALTITUDE  DEPENDENCE  OF  THE  FREEZING  SCALE  COMPARISON  OF 
MEASUREMENTS  WITH  EXISTING  THEORY 


As  was  discussed  in  Section  4,  to  make  the 
upleg  and  downleg  PLUMEX  data  consis¬ 
tent  in  terms  of  a  unique  altitude 
dependence  of  the  freezing  scales,  it  was 
necessary  to  map  the  data  to  equal  altitudes 
at  the  magnetic  equator.  This  seems 
reasonable  because  that  is  where  the  driver 
is  concentrated.  We  next  wish  to  consider 
how  the  altitude  dependence  found  com¬ 
pares  with  existing  theoretical  models. 

One  prominent  model  for  the  freezing 
scale,  due  to  Sperling  and  Glassman  (1985) 
has  a  relatively  simple  altitude  dependence. 
In  that  theory,  the  freezing  scale  is 
controlled  by  a  balance  between  growth  and 
decay  mechanisms,  and  has  a  uirr'  2 
dependence  with  altitude,  where  u,„  is  the 
ion-neutral  collision  frequency.  The  Sperl¬ 
ing  and  Glassman  (1985)  prediction,  for 
exospheric  temperature  of  800  K  is 
superimposed  in  Figure  18.  The  agreement 
between  observation  and  theory  is  reason¬ 
able.  particularly  at  the  lowest  and  highest 
altitudes.  Overall,  however,  the  slope  of  the 
prediction  curve  is  too  steep  to  match  the 
bulk  of  the  points.  This  suggests  that  a 
linear  dependence  on  Ujn“ *  might  provide  a 
better  fit  to  the  break  scale  data,  at  least 
for  the  nonbottomside  samples.  This  is,  in 
fact,  the  case,  as  the  curve  (labeled 
magnetic  equator)  in  Figure  18  illustrates. 
Hie  linear  growth  rate  is  equal  to  y  = 
gcn/i>i„L.  where  g fff  is  the  effective 
gravitational  constant  and  L  is  the  gradient 
scale  length.  Using  a  simple  dipole 
magnetic  field  model  and  taking  into 
account  the  angle  that  the  gravitational 


vector  makes  with  the  field,  it  can  be  shown 
that  gefr  is  nearly  constant  along  a  given 
field  line.  The  observed  relationship 
between  the  break  wavelength  and  uin 
altitude  profiles  suggests  a  dependence  of 
the  break  wavelength  on  -ycff  =  gcff/uin. 
independent  of  the  gradient  scale. 

To  accept  a  Vj,,'1  dependence  of  the  spectral 
break  over  the  central  portion  of  the  F 
layer,  some  explanation  is  necessary  as  to 
why  the  bottomside  spectral  forms  do  not 
fit  the  same  pattern.  All  of  the  spectra  that 
deviate  from  the  vin  curve  at  low  altitudes 
occur  below  about  280  km  true  altitude 
and,  we  note  from  Figure  13,  these 
segments  are  all  on  the  F-  layer  bottomside. 
In  terms  of  the  overall  data  set.  these  points 
are  different  from  those  at  higher  altitudes: 
each  is  on  a  steep  positive  gradient  and 
exhibits  a  large  value  of  8N/N.  This  is 
illustrated  more  specifically  in  Figure  19. 
which  is  a  scatter  plot  of  gradient  versus 
oN/N.  with  the  bottomside  points  circled. 
Both  the  steep  gradient  and  the  large  oN/N 
are  indicators  that,  although  the  upper  F 
layer  may  be  stable,  the  bottomside  may  be 
actively  structuring.  If  local  conditions 
(rather  than  those  at  the  magnetic  equator) 
dominate,  the  theoretical  curves  in  Figure 
19  may  move  downward  as  much  as  60  km 
in  altitude.  This  local  structuring  vin' 1 
dependence  is  also  shown  in  Figure  19.  and 
provides  a  reasonably  good  fit  to  the 
bottomside  dependence  we  observe.  This 
may  be  a  fortuitous  agreement,  however, 
and  it  should  be  kept  in  mind  that  several 
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GRADIENT  SCALE 


other  factors  could  explain  the  difference  possibility  of  a  vertical  velocity  shear 

in  break  wavelength  behavior  on  the  ( Tsunoda  et  at..  1981),  and  the  latitudinal 

bottomside.  Included  among  these  is  the  dependence  of  the  eastward  neutral  wind. 


NORMALIZED  STND  DEV  DENSITY 


Figure  iy.  Scatter  Plot  of  Gradient  Scale  versus  oN/N  measured  by  PLUMEX. 
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SECTION  6 


CONCLUSIONS 


In  this  report,  we  have  made  a  comprehen¬ 
sive  experimental  study  of  the  spectral 
characteristics  of  equatorial  ionospheric 
F-region  structure  Included  are  the  charac¬ 
terization  of  spectra  measured  using  RF 
propagation,  and  in  situ  satellite  and  rocket 
observations  scanning  both  vertically  and 
horizontally  across  the  magnetic  field. 

The  vertical-scan  PLUMEX  rocket  data  are 
the  easiest  to  interpret  in  terms  of  the 
physical  processes  that  produce  the  irregu¬ 
larities.  This  is  because  the  structuring 
processes  appear  to  be  consistent  oxer  the 
time  and  spatial  extent  of  the  measure¬ 
ment.  Under  those  conditions,  the  irregu¬ 
larity  power  spectra  on  the  topside  of  the 
F  layer  (i.e..  inertia-dominated  altitudes) 
are  power  law  in  form,  with  a  one-dimen¬ 
sional  slope  p  -  2.  Lower  in  altitude  (i.e.. 
collision  dominated),  the  spectra  are  two 
component  in  form,  with  a  shallow 
intermediate  scale  regime  and  a  more 
steeply  sloped  dissipation  regime.  The 
freezing  scale,  which  we  have  associated 
with  the  break  of  slope  in  the  two-compo¬ 
nent  spectra,  moves  to  larger  scale  sizes 
with  increasing  altitude.  This  altitude 
dependence  is  well  matched  by  the  inverse 
of  the  ion-neutral  collision  frequency.  The 
comparison  of  the  upleg  and  downleg  data 
suggests  that  the  freezing  scale  characteris¬ 
tics  of  an  irregularity  volume  is  determined 
bv  local  conditions  at  the  time  and  location 
of  their  origin. 

The  AE-E  satellite  data  that  have  been 
analyzed  were  collected  over  a  large  spatial 
volume  and  over  many  hours  of  local  time. 
Approximately  one-third  of  those  data  are 
single  power  law  in  form,  which  is 


consistent  with  the  satellite  altitude  and  the 
PLUMEX  observations.  Also  consistent 
with  the  PLUMEX  data  is  an  average  single 
power-law  slope  of  p  ~  2.  Among  the  domi¬ 
nant  two-component  population,  there  is  no 
obv  ious  collision  dependence  of  the  freez¬ 
ing  break,  even  with  changing  satellite 
altitude.  In  light  of  the  PLUMEX  data  and 
in  particular,  the  implication  that  the 
spectral  shape  is  not  necessary  a  reflection 
of  local  conditions,  this  is  not  surprising. 
On  the  other  hand,  the  average  PLUMEX 
vertical  and  AE-E  horizontal  spectral 
shapes  are  in  good  agreement:  this  suggests 
that  spectral  anisotropy,  if  present,  is  not 
major.  The  average  break  wavelengths  of 
the  vertical  spectra  are  somewhat  larger 
than  the  AE-E  average,  but  are  well  w  ithin 
the  statistical  distribution  observed  by  the 
satellite.  The  shallow  intermediate  slopes 
seen  by  AE-E  anti  PLUMEX  are  essentially 
the  same  with  a  value  p;i  ~  0.8.  The  dissi¬ 
pation  regime  average  slopes  from  AE-E 
show  a  significant  change  with  altitude,  and 
are  shallower  than  those  seen  by  PLUMEX 
or  implied  from  tf^phase  scintillation  data. 
The  current  analysis  has  not  addressed  this 
difference  in  detail,  but  subtle  spectral 
anisotropy  seems  likely  to  be  the  cause.  A 
value  of  pi,  =  3.2  can  be  used  as  an  esti¬ 
mated  average  overall  value  of  vertical/ 
horizontal,  one-dimensional  slope  in  the 
dissipation  regime. 

Care  must  be  taken  to  avoid  strong  scatter 
effects  if  phase  scintillation  data  are  to  be 
used  to  imply  the  spectral  form  of  in  situ 
irregularities.  In  addition,  a  significant  per¬ 
centage  of  the  phase  spectra  appear  to  be 
altered  by  superposition  of  energy  from 
different  altitudes.  This  is  inherent  in  the 
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technique,  but  can  be  avoided  by  careful 
analysis.  These  disadvantages  of  the 
method  are  more  than  offset  by  its  ease  and 
economics  as  a  morphological  tool.  The 
slopes  of  the  single  power-law  phase  data 
agree  very  well  with  the  slopes  observed  in 
situ.  The  average  value  is  near  p  =  3,  which 
translates  to  a  onedimensional  slope  of  p 
=  2.  This  suggests  that  thin  phase  screen 
theory,  which  predicts  this  correspondence 
in  weak  scatter,  is  accurate.  The  two- 


component  phase  spectral  form  is.  on 
average,  steeper  than  the  in  situ  spectra  in 
the  intermediate-scale  regime.  The  varia¬ 
tion  of  this  slope  between  station  longitudes 
suggests  that  it  is  a  result  of  leakage  of 
superposed  phase  energy  from  structures 
much  larger  than  a  few  kilometers.  The 
dissipation  regime  slopes  are,  as  noted 
above,  in  excellent  agreement  with  those 
seen  by  PLUMEX,  but  steeper  than  the 
average  in  the  AE-E  statistics. 
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HARRY  DIAMOND  LABORATORIES 

ATTN:  SLCIS  IM-TL  (TECH  LIB) 


INFORMATION  SYSTEMS  COMMAND 

ATTN:  STEWS-TE-N  K  CUMMINGS 

U  S  ARMY  ATMOSPHERIC  SCIENCES  LAB 
ATTN:  SLCAS-AE  (DR  NILES) 

ATTN:  SLCAS-AE-E 
ATTN:  SLCAS-AR  DR  H  HOLT 

U  S  ARMY  FOREIGN  SCIENCE  &  TECH  CTR 
ATTN:  DRXST-SD 

U  S  ARMY  MISSILE  COMMAND 

ATTN:  AIAMS-S/B  J  GAMBLE 

U  S  ARMY  MISSILE  COMMAND/AMSMI-RD-CS-R 
ATTN:  AMSMI-RD-CS-R  (DOCS) 

U  S  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
ATTN:  MONA-NU 

U  S  ARMY  NUCLEAR  EFFECTS  LABORATORY 
ATTN:  ATAA-PL 

U  S  ARMY  STRATEGIC  DEFENSE  CMD 
ATTN:  CSSD-H-LS  B  CARRUTH 
ATTN:  CSSD-H-SA 
ATTN:  CSSD-H-SA/R  SMITH 
ATTN:  CSSD-H-SAV 
ATTN:  CSSD-IN-T  M  POPE 

U  S  ARMY  STRATEGIC  DEFENSE  COMMAND 
ATTN:  ATC-OW  DAVIES 
ATTN:  CSSD-GR-SW  DICKINSON 

USA  SURVIVABILITY  MANAGMENT  OFFICE 
ATTN:  SLCSM-SE  J  BRAND 

DEPARTMENT  OF  THE  NAVY 

NAVAL  ELECTRONICS  ENGRG  ACTVY,  PACIFIC 
ATTN:  CODE  250  DOBRYHIM 

NAVAL  OCEAN  SYSTEMS  CENTER 

ATTN:  CODE  542,  J  FERGUSON 

NAVAL  RESEARCH  LABORATORY 

ATTN:  CODE  2000  J  BROWN 
ATTN:  CODE  2627  (TECH  LIB) 

2  CYS  ATTN:  CODE  4100  H  GURSKY 

ATTN:  CODE  4121.8  H  HECKATHORN 

ATTN:  CODE  4183 

ATTN:  CODE  4700  S  OSSAKOW 

ATTN:  CODE  4701 

ATTN:  CODE  4720  J  DAVIS 

ATTN:  CODE  4780  B  RIPIN 

ATTN:  CODE  4780  DR  P  BERNHARDT 

ATTN:  CODE  4780  J  HUBA 

ATTN:  CODE  5300 

ATTN:  CODE  5326  G  A  ANDREWS 

ATTN:  CODE  5340  E  MOKOLE 

ATTN:  CODE  8344  M  KAPLAN 

NAVAL  TECHNICAL  INTELLIGENCE  CTR 
ATTN:  DA44 
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NAVAL  UNDERWATER  SYSTEMS  CENTER 
ATTN:  CODE  3411,  J  KATAN 

OFFICE  OF  CHIEF  OF  NAVAL  OPERATIONS 
ATTN:  OP941D 

OFFICE  OF  NAVAL  RESEARCH 

ATTN:  1132SM/A  TUCKER 

SPACE  &  NAVAL  WARFARE  SYSTEMS  CMD 
ATTN:  PD50TD 
ATTN:  PD50TD1  G  BRUNHART 
ATTN:  PME  106-4  S  KEARNEY 
ATTN:  PME-106  F  W  DIEDERICH 

THEATER  NUCLEAR  WARFARE  PROGRAM  OFC 
ATTN:  PMS-42331F  (D  SMITH) 

DEPARTMENT  OF  THE  AIR  FORCE 

AFIA/INIS 

ATTN:  AFIA/INKD  MAJ  SCHROCK 

AIR  FORCE  CTR  FOR  STUDIES  &  ANALYSIS 
ATTN:  AFCSA/SAMI 
ATTN:  AFCSA/SASC 

AIR  FORCE  GEOPHYSICS  LABORATORY 
ATTN:  J  KLOUBACHAR 
ATTN:  OP/WBLUMBERG 
ATTN:  SANTI  BASU 
ATTN:  SUL 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RSCH 
ATTN:  AFOSR/NP 

AIR  FORCE  SPACE  SYSTEMS  DIVISION 
ATTN:  YA 
2CYS  ATTN:  YN 

AIR  FORCE  TECHNICAL  APPLICATIONS  CTR 
ATTN:  TN 

AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL-LSE 

STRATEGIC  AIR  COMMAND/XRFS 
ATTN:  XRFS 

WEAPONS  LABORATORY 
ATTN:  NTCA 
ATTN:  WL/SUL 

DEPARTMENT  OF  ENERGY 

EGAG,  INC 

ATTN:  D  WRIGHT 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  L-97  T  DONICH 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  RW  WHITAKER 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  A  D  THORNBROUGH 
ATTN:  CODE  9014  R  BACKSTROM 
ATTN:  D  DAHLGREN 
ATTN:  DIV  2344  ROBERT  M  AXLINE 


ATTN:  ORG  1231  JR  LEE 
ATTN:  ORG  91 14  WD  BROWN 
ATTN:  TECH  LIB 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  OSWR/NED 
ATTN:  OSWR/SSD  FOR  L  BERG 

DEPARTMENT  OF  COMMERCE 
ATTN:  G  REEVE 
ATTN:  J  HOFFMEYER 
ATTN:  WUTLAUT 

U  S  DEPARTMENT  OF  STATE 
ATTN:  PM/TMP 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN:  A  MORSE 
ATTN:  BRIAN  PURCELL 
ATTN:  C  CREWS 
ATTN:  CRICE 
ATTN:  G  LIGHT 
ATTN:  IGARFUNKEL 
ATTN:  J  KLUCK 
ATTN:  MROLENZ 

ATMOSPHERIC  AND  ENVIRONMENTAL  RESEARCH  INC 
ATTN:  MKO 

AUSTIN  RESEARCH  ASSOCIATES 
ATTN:  J  THOMPSON 

BERKELEY  RSCH  ASSOCIATES,  INC 
ATTN:  CPRETTIE 
ATTN:  J  WORKMAN 
ATTN:  NTGLADD 
ATTN:  S  BRECHT 

CALIFORNIA  RESEARCH  A  TECHNOLOGY,  INC 
ATTN:  M  ROSENBLATT 

CHARLES  STARK  DRAPER  LAB,  INC 
ATTN:  ATETEWSKI 

CORNELL  UNIVERSITY 

ATTN:  D  FARLEY  JR 
ATTN:  M  KELLY 

DEFENSE  GROUP,  INC 

ATTN:  MK  GROVER 

ELECTROSPACE  SYSTEMS,  INC 
ATTN:  P  PHILLIPS 

EOS  TECHNOLOGIES,  INC 
ATTN:  B  GABBARD 
ATTN:  R  LELEVIER 

GENERAL  ELECTRIC  CO 

ATTN:  ROBERT  H  EDSALL 

GENERAL  RESEARCH  CORP  INC 
ATTN:  JEOLL 
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GEO  CENTERS,  INC 

ATTN:  E  MARRAM 

HARRIS  CORPORATION 
ATTN:  EKNICK 

HSS,  INC 

ATTN:  D HANSEN 

INFORMATION  SCIENCE,  INC 
ATTN:  WDUDZIAK 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  E  BAUER 

J  S  LEE  ASSOCIATES  INC 
ATTN:  DR  J  LEE 

JAYCOR 

ATTN:  J  SPERLING 

JOhNS  HOPKINS  UNIVERSITY 
ATTN:  CMENG 

KAMAN  SCIENCES  CORP 
ATTN:  DASIAC 
ATTN:  E CONRAD 
ATTN:  G  DITTBERNER 

KAMAN  SCIENCES  CORPORATION 
ATTN:  BGAMBILL 
ATTN:  DASIAC 
ATTN:  R RUTHERFORD 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  J  HENLEY 
ATTN:  J  KUMER 
ATTN:  R SEARS 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN.  D  T  RAMPTON 

LTV  AEROSPACE  &  DEFENSE  COMPANY 
2  CYS  ATTN:  LIBRARY 

M  I  T  LINCOLN  LAB 

ATTN:  D  TOWLE 
ATTN:  M  LEE 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN:  J  BENNETT 

MAXIM  TECHNOLOGIES,  INC 
ATTN:  J  SCHLOBOHM 

MCDONNELL  DOUGLAS  CORP 
ATTN:  TCRANOR 


MISSION  RESEARCH  CORPORATION 


ATTN: 

J  KENNEALY 

ATTN: 

R  ARMSTRONG 

ATTN: 

S  DOWNER 

ATTN: 

W  WHITE 

MISSION  RESEARCH  CORP 

ATTN: 

B  R  MILNER 

ATTN: 

C  LONGMIRE 

ATTN: 

D ARCHER 

ATTN: 

D  KNEPP 

ATTN: 

D  LANDMAN 

ATTN: 

DOCUMENT  CONTROL/R  BOGUSCH 

ATTN: 

F  FAJEN 

ATTN: 

F  GUIGLIANO 

ATTN: 

G  MCCARTOR 

ATTN: 

K  COSNER 

ATTN: 

M  FIRESTONE 

ATTN: 

R  BIGONI 

ATTN: 

R  DANA 

ATTN: 

R  HENDRICK 

ATTN: 

RKILB 

ATTN: 

S  GUTSCHE 

ATTN: 

TECH  LIBRARY 

MITRE  CORPORATION 

ATTN:  DRMRDRESP 

MITRE  CORPORATION 

ATTN:  RCPESCI 
ATTN:  W  FOSTER 

NORTHWEST  RESEARCH  ASSOC,  INC 
ATTN:  E  FREMOUW 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  HBRODE 

PHOTOMETRICS,  INC 

ATTN:  ILKOFSKY 

PHOTON  RESEARCH  ASSOCIATES 
ATTN:  D  BURWELL 
ATTN:  0  LEWIS 

PHYSICAL  RESEARCH  INC 
ATTN:  WSHIH 

PHYSICAL  RESEARCH  INC 
ATTN:  ACECERE 

PHYSICAL  RESEARCH  INC 
ATTN.  HFITZ 
ATTN:  PLUNN 
ATTN:  WHEUSER 


MCDONNELL  DOUGLAS  CORP 
ATTN:  J  GROSSMAN 
ATTN:  RHALPRIN 

METATECH  CORPORATION 
ATTN:  R  SCHAEFER 
ATTN:  WRADASKY 

METEOR  COMMUNICATIONS  CORP 
ATTN:  R LEADER 


PHYSICAL  RESEARCH,  INC 
ATTN:  R  DELIBERIS 
ATTN:  T  STEPHENS 

PHYSICAL  RESEARCH,  INC 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  W  SCHLUETER 

PHYSICS  INTERNATIONAL  CO 
ATTN:  C  GILMAN 
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R  &  D  ASSOCIATES 

ATTN:  F  GILMORE 
ATTN:  MGANTSWEG 

R  &  D  ASSOCIATES 

ATTN:  J  WALTON 

RANDCORP 

ATTN:  C  CRAIN 
ATTN:  E  BEDROSIAN 

RAND  CORP 

ATTN:  B BENNETT 

RJO  ENTERPRISES/POET  FAC 

ATTN:  STEVEN  KRAMER 
ATTN:  W BURNS 

S-CUBED 

ATTN:  C NEEDHAM 
ATTN:  T CARNEY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  S  ROSENCWEIG 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  C  SMITH 
ATTN:  D  HAMLIN 
ATTN:  D SACHS 
ATTN:  LLINSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J COCKAYNE 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  HSUNKENBERG 
ATTN:  LIBRARY 


SRI  INTERNATIONAL 
2  CYS  ATTN:  J  F  VICKREY 
2  CYS  ATTN:  R  LIVINGSTON 
ATTN:  R  T  TSUNODA 
ATTN:  W  CHESNUT 
ATTN:  WJAYE 

STEWART  RADIANCE  LABORATORY 
ATTN:  R  HUPPI 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  J  WOLFSBERGER,  JR 
ATTN:  N  PASSINO 

TOYON  RESEARCH  CORP 
ATTN:  J  ISE 

TRW  INC 

ATTN:  DRDGRYBOS 

TRW  SPACE  &  DEFENSE  SYSTEMS 
ATTN:  DM  LAYTON 

USER  SYSTEMS,  INC 

ATTN:  S  W  MCCANDLESS,  JR 

VISIDYNE,  INC 

ATTN:  J CARPENTER 

DIRECTORY  OF  OTHER 

BOSTON  UNIVERSITY 

ATTN:  MICHAEL  MENDILLO 
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